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1. Introduction 



The one-dimensional Schrodinger equation 

Lip(x) = —ip"(x) + u(x)ip(x) = Xip(x) (1) 

plays an important role in the theory of nonlinear evolution equations, where it serves as an 
auxiliary spectral problem allowing to integrate the KdV-equation. From this point of view 
the KdV-evolution of u(x) in time yields the isospectral deformations of ([!]). In the quantum 
mechanical context, solutions of (|l|) subjected to some boundary conditions describe physical 
states of a particle moving on the line. In both applications a special attention is drawn to the 
class of exactly solvable potentials u(x) whose discrete spectra are given by some known functions 
of a quantum number. The powerful approach to derivation of such potentials is provided by the 
factorization method |]J based on the Darboux transformations. In this method one has a set of 
operators Lj = —d 2 /dx 2 + Uj(x) which are factorized as products of linear differential operators 
Alj = —d/dx + fj(x), Aj = d/dx + fj(x) up to some constants A.,-: Lj = A^ Aj + Aj. Since the 
spectra of the operators A^Aj and Aj Aj may differ only by one lowest level (provided fj has 
no severe singularities), the Hamiltonians Lj + i and Lj will have close spectral properties if we set 
Lj + i = A~j +1 Aj +1 + Xj+i = AjAj + Aj. This constraint is equivalent to the following differential 
equation relating potentials Uj and 

/;(*) + fj +1 (x) + fji.r) - ff +1 (x) = Xj . , - A/. (2) 

The chain of such equations is called the dressing chain in the theory of solitons. It is quite useful 
for studying symmetries of the Schrodinger equation. For particular solutions of the dressing chain 
parameters Aj coincide with ordered discrete spectra of the corresponding Hamiltonians Lj. So, 
the spectra of harmonic oscillator, Coulomb, and other "old" solvable potentials are easily found 
after the plugging in @ the ansatz: fj(x) = a(x)j + b(x) + c(x)/j 0. 

In this Letter we discuss the particular class of the self-similar potentials related to each other 
by the scaling of argument: Uj+\{x) = q 2 Uj(qx), which emerges as a result of the constraints 
fj + i(x) = qfj(qx), Xj + i = q 2 Xj imposed on (0). This class has been studied by A.Shabat for 
the real x and q, < q < 1, when it corresponds to "infmite-soliton" systems, i.e. when the 
potentials are reflectionless and decrease slowly at space infinities. The corresponding discrete 
spectra are purely exponential; they are generated by the g-deformed Heisenberg-Weyl, or q- 
oscillator algebra 0. In the two limiting cases, q — >• and q — > 1, such potentials are reduced to 
the one-soliton and harmonic oscillator potentials respectively. In the present paper we consider 
this self-similar system in the complex domain of the coordinate x and parameter q. First we 
analyze the existence and uniqueness of the corresponding functions fj(x). Then we characterize 
the qualitative structure of singularities exhibited by fj(x). When q is a primitive n-th root of unity 
the situation is simplified. The odd n cases are shown to be related to the finite-gap potentials. 
The even n cases may contain a functional non-uniqueness. The general q 3 = 1 solution and 
a special q 4 = 1 solution are expressed through the equianharmonic and (pseudo)lemniscatic 
Weierstrass functions respectively. Algebraically, all these root-of-unity potentials are naturally 
related to the representations of the g-oscillator algebra at q n = 1. A more wide class of the 
self-similar potentials leading to more complicated g-deformed algebras has been described in f|. 
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Though some of our results may be reformulated easily for the latter systems, they will not be 
discussed here. 

Appearance of the g-analogs of the harmonic oscillator and other spectrum generating alge- 
bras raises special interest to the self-similar potentials. This is inspired by the recent inten- 
sive discussion of quantum algebras, or g-deformations of Lie algebras ||, which bear universal 
character due to the large number of applications. The g-oscillators themselves were reinvented 
as the related objects ||. It was found that a natural group-theoretical setting for the basic, 
or g-hypergeometric functions |7[] and corresponding orthogonal polynomials is provided by the 
quantum algebras when they are realized with the help of purely finite-difference operators j5|. 
However, the g-hypergeometric functions are defined in general only at \q\ < 1 and we were not 
able to find in the literature an example of g-special function for q n = 1 - the cases when the rep- 
resentation theory of quantum algebras essentially differs from that for the standard Lie algebras 
IJ. The self-similar potentials provide a realization of the creation and annihilation operators 
of the g-oscillator algebra in terms of both differential and finite-difference operators. It is the 
presence of the differential part that principally differs our approach from the mentioned above 
and leads to the well-defined g-special functions at q n = 1. 

The paper is organized as follows. In Sect. 2 we describe the needed facts about the periodic 
dressing chain. In Sect. 3 we present the class of self-similar potentials and discuss its analytical 
properties. The central result here is the theorem on the existence and uniqueness of solutions for 
the basic differential equation with deviating argument (|14|) when |g| < 1, or q n = 1. The g 3 = 1 
and g 4 = 1 cases are considered in detail in Sects. 4 and 5 respectively. Sect. 6 is devoted to 
the description of g-oscillator algebra representations related to the taken g-transcendent. In the 
Appendix we present exact solution of the spectral problem associated with the g 3 = 1 system. 

The results obtained in this paper were partially presented at the Canadian Mathematical 
Society Meeting (Montreal, December 1992). 

2. The Periodic Dressing Chain 

We discuss here some results of the theory of nonlinear evolution equations on the basis of the 
dressing chain. 

Following PD |, one can rewrite (|1]) as a matrix equation and introduce a chain of associated 



spectral problems 
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V j (x) = U j * j (x), U j =[^_ x Q j, *i=(^J> 3^Z ( 3 ) 
solutions of which are related to each other by the transformations 

V :(•'•) = ''/I', (•'')• (4) 

where Bj is a 2 x 2 matrix whose entries are polynomials of the spectral parameter A. The 
compatibility condition for (|)-(|) leads to the equation 

B> = I ) . , n, - B 3 U y (5) 
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When the map ipj — > 4>j+i does not depend on A the solution of fl5|) looks as follows 



B j 



fj{x) 1 
fji-r) • A; • A fj ( x ) 



(6) 



U A X ) = fj - fj + A i, u j+ x-Uj = 2f-, (7) 

where A.,- are some constants of integration. Substituting the first of relations (|7|) into the second 
one we get the dressing chain (0). 

Denote by A, = 5j + jy-i • • • Bj + xBj the product of Darboux transformation matrices (^|). From 
(|Sp it follows that 

4 = U j+N Aj - AjUj. (8) 

This equation allows to integrate the chain (0) for the special class of potentials defined by the 
periodicity condition 

U j+N = U j , or f j+N (x) = fj(x), A j+A r = A j , (9) 

which describes a finite-dimensional dynamical system. For the odd N this leads to the finite-gap 
potentials flOfl . 

For convenience we fix the index j in fl8|), (|j), j = 0, and set L = L , U = Uq, A = A . By 
its definition the matrix A maps solutions of the corresponding Schrodinger equation onto each 
other and, hence, satisfies the matrix Lax equation A = [U,A]. From the latter, one can derive 
the following relation for the matrix element a\2 = 

|/3/3" + (A - u )f3 2 - J/3' 2 = det A - \{Tr Af . (10) 

Both Tr A and det A are easily seen to be integrals of motion; at iV = 2k + 1 one has 

2k 

detA = n(A-Ai) , Tr A = (-l) k (I \ k + h\ k - 1 + ■ ■ ■ + I k ). 
Equation (|T0|) can be rewritten now in the form 

u -X = F 2 -F\ F(x) = ^^, (11) 

where 

2k 



W 
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{TrAf - AdetA = -4 JJ(A - EA. 



i=0 

If all the constants E\ are real and mutually different then the spectrum of non-singular potentials 
Uj(x) has zonal structure and Ei coincide with the boundaries of gaps. It is well-known that such 



potentials can be expressed in terms of the Riemann G-function [11 



d 2 

u(x) = -2— In e(l(x)) + const, (12) 
dx 2 

k k 

0(0 = ex p{ H m p b ps m s + l P m p)i 

m 1 ,...,m k £Z P,s=l p=l 
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where l p (x) are linear functions of x whose coefficients together with the parameters b ps are de- 
termined by Ei. 

We shall not go further in the presentation of general results. For a detailed account of 
integrability of the dressing chain at even period N and in the case of a more general than @ 
closure condition Uj+N = Uj + a, where a is a constant, we refer to the recent work O . 



3. Self-Similar Potentials 

The particular solution for the dynamical system (0), which we investigate in the present paper, 
is fixed by the following self-similarity condition 

fj( x ) = tff (fx), fij = g 2 >, Hj = Aj+i - Xj. (13) 

In this case an infinite number of relations (@) is reduced to one differential equation with the 
deviating argument which was introduced by A.Shabat in f\ : 

Y x (f( x ) + Qffa)) + f 2 W ~ 9 2 / 2 (gx) = /i. (14) 



Consider the initial value problem for the equation ([14]) : 

f(0) = a < oo. (15) 

We fix the initial condition at the point x = because it is a fixed point of scaling x — > qx (there 
is another such point x = oo which was considered in ||14|| ). The parameter /x is redundant - it 
may be removed by rescaling of the coordinate and f(x), but we prefer to keep it as a unique 
dimensional parameter. 

A more general class of the self-similar potentials arises after the following g-periodic closure 
of the dressing chain Q : 

f j+N (x) = qfj(qx), fi j+N = q 2 fij. (16) 

Here we limit our discussion to the N = 1 case which coincides with (|l3|) . 

Let f(z, q) denotes a function of two complex variables that satisfies the equation (|i~4D. Plug- 
ging into (|i~4|) the power series expansion 

oo 

/(z,g) = $>„(g)z", (17) 

n=0 

one obtains the recursion formula for the coefficients cL: 



n+2 _ i 

~,n+2\~ _ y 1 



i + q n+2 ya n+ i = H —— E 

ai(l + g 2 ) =fi+ (q 2 - l)o 



n > 1, 



s=0 

2\ i / 2 i \~2 



3 According to the classification of equations with a deviating argument |13|, equation (1l4) is of the neutral type. 
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When f(0,q) = 0, all with even k vanish, i.e. the f{z,q) is an odd function of z which 
corresponds to symmetric potentials.^ For this case it is convenient to rewrite formulae (|T7|) , (|18|) 
as follows: 

oo 

f(z,q) = J2an(q)z 2n -\ (19) 

n=l 

(1 + q 2n )a n = V a s a n - s , a Y = — — -. (20) 

2n — 1 ^ 1 + <r 

Lemma. Series ([HJ), (|20|) converges in the disc \z\ < R q , for every fixed value of q, \q\ < 1. 



For the radius of convergence R n the following estimate holds: R„ > — Z — , where a = ] + \ q \ 2 ■ 

to q to 9 - 2 v /|a 1 |a' !-kl 

Proof. For any natural number p and |g| < 1 one has: | ^ q i v \ < a. As a result, \a n \ < c ra , where 
c n -coefficients are defined by the same recursion relations (^) with ci = |ai| and the ? 2 "^} -factor 
being replaced by the a for all n. But the series J2™=i c n \z\ 2n ~ l is proportional to a scaled tangent 
function of \z\. Hence, our series is majorized by 



\f(z,q)\ < ^tanyVil" \ z \ 

which asserts the statement of the lemma. At q = this gives an exact answer for R q . Using 
analogous method one can prove convergence of the series ([17]) near z = as well. □ 



Solutions of (|14D for |g| > 1 can be obtained from those for \q\ < 1 by the following relation: 
f(z,q~ x ) = iqf(—tqz,q); therefore it is sufficient to consider q from the unit disk: \q\ < 1. 

Theorem. For any q such that \q\ < 1, or q 2k+1 = 1, k — 0, 1, 2, ... in some neighbourhood of 
z = there exists the unique solution f(z,q) of ( |H| ) satisfying the initial condition (13). When 



q is a primitive root of unity of even degree, q 2k = 1, existence of analytic solutions depends on 
the value of ao; depending on k the solution satisfying properly chosen initial condition may be 
non-unique. 

Proof. Consider separately the following cases: 

1. Suppose that \q\ < 1. Then, by the given f(0,q) one can determine uniquely f'(0,q), 
f"(0,q), etc by taking successive derivatives of fll4]) . Therefore the formal series ( [P7D is 
defined uniquely, and by the lemma it converges near the z = point. So, solution exists 
and it is unique.Q 

2. Let q be a primitive root of unity of odd degree: q 2k+1 = 1. In this case recursion relations 
(|18D, (|20| ) define uniquely series coefficients. Their convergence in some neighborhood of 
z = follows from the lemma because the lower bound for R q can be estimated by replacing 
a — > at, where is defined by the inequality < cik for any natural p and q = exp T^rpi- 
Since q n ^ — 1 such finite does exist but its value depends on k. 



6 At least for |g| < 1, as it will be seen. 

7 For real q, < q < 1, this statement follows from a general theorem on existence and uniqueness of solutions 
of differential equations with deviating argument of the neutral type proved by G.A. Kamenskii |13|. 



5 



3. Let q be a primitive root of unity of even degree: q 2k = 1. Now the left hand side of the 



recursion formula fll8|) vanishes at n + 2 = fc(2Z + 1), Z = 0, 1, . . ., but the right hand side 
does so only for special values of the parameter clq. For these specific initial conditions 
the coefficients dk(2i+i)-i of the series (|17| ) can be defined in a non-unique way (i.e. some, 
or all of them may take arbitrary values), any other choice of So leads to the mismatch 
in recursion relations for a n and therefore (|T4]) has no solutions analytical at z — 0. For 
example, when q = —1, fi ^ 0, the analytic solution is unique and exists only for do = 
(it is f(z, —1) = fiz/2). If = then there is no restriction on a and any even function 
f(z) = f(—z) is a solution; the recursion formula ( |I%D gives in this case an ambiguity in each 
a,2n- When q 2 = —1, an analytic solution exists iff 2d\ = //. Then, the derivative f'(Q,q) 
is not determined by the equation flliD and it can take any value. In particular, one has 
f(z, ±z)a =^=o = d\z + a 5 z 5 + dgz 9 + . . . , where d\, S5, . . . are arbitrary coefficients such that 
the series converges. In section 5 below we give an exact dependence of the solutions on an 
arbitrary function for this case. Similarly, when g 3 = —1 one should put d = 0, ^ r,J 



, u — 

and ambiguous coefficients are 0,2, S 8 , du, .... □ 

Let us make two remarks. First, for the general system of equations (|2J), (0) with N = 2n the 
uniqueness of its solutions may be violated if Iq = 2 J2iLi fi — ( see 10)- This condition is obvi- 
ously satisfied in the self-similar case when q N = 1, but there are also other situations exhibiting 
non-uniqueness. Second, analysis of the existence of analytical solutions f(z, q) when q = e 2mr ^, 
- an irrational number, is essentially a number theory problem. The nature of irrationality of 
plays a crucial role in the determination of the asymptotic growth of the series coefficients, in 
this respect our problem is close to the question on the convergence of g-hypergeometric series at 
such q. 

Since at \q\ < 1 the function f(z,q) is analytic in the neighbourhood of z = 0, one can in 
principle construct f(z, q) on the whole complex plane as follows. Let us rewrite equation (|TJ|) 
in the form f'(z) + f 2 (z) = v(qz), where v(z) = q 2 (f 2 {z) — f'{z)) + fi, and fix some number 
R, < R < R q . In the coordinate region R < \z\ < R/\q\ this is just the Riccati equation 
because v(qz) is the fixed function determined by the values f(z,q) at \z\ < R. The general 
solution of this equation is: f(z,q) = d/dzln (w\ + cw 2 ), where u>i )2 are independent solutions 
of the auxiliary Schrodinger equation —w"(z) + v(qz)w(z) = 0. The constant c is fixed by some 
boundary condition. Continuing this procedure iteratively in the rings R/\q\ n < \z\ < -R/|g| n+1 , 
we recover f(z, q) on the open complex z-plane. Since all singularities of f(z, q) appear from zeros 
and singularities of the solutions of a linear differential equation, one can find their qualitative 
structure. 

Proposition 1. The function f(z, q) is meromorphic in the open complex z-plane for any fixed 
q, \q\ < 1. 

Proof. Because v(z) is analytic in the disc \z\ < R, all permitted singularities of f(z, q) in the 
first ring R < \z\ < R/\q\ are simple poles with the residues r\ — 1. These poles originate from 
simple zeros, w(z p ) = 0, of the wave function w(z) fixed by the boundary condition w'/w\ z= r = cq. 
In the ring R < \z\ < R/\q\ the function v(z) may contain double poles at the points z = z p /q 
with the residues equal to ri(ri + 1) = 2. There may be also accompanying simple poles, but they 
are irrelevant. In the vicinity of these singularities solutions of the auxiliary Schrodinger equation 
have the form w(z) oc (z — z p /q) r where r is found from the indicial equation r(r — 1) = riiji + 1). 
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This equation has two solutions: r = — n and r = r% + 1 corresponding to linearly independent 
functions ifi^- The boundary condition f(R/\q\,q) = C\ fixes their relevant combination. Since 
u>i.2 may have only singularities (or zeros) of the form (z — z p /q) r and simple zeros, we conclude 
that the function f(z,q) in the second ring can have only pole- type singularities with the integer 
residues r% — r, or 1. Continuing this procedure iteratively from one ring to another we prove the 
assertion. It is convenient to call the poles with the residues equal to +1, which arise due to the 
simple zeros of w(z) and are located at z = z p , as the primary ones. The poles arising from them 
as the descendents in the points z = z p /q n with the residues equal to one of the possible values 
of r G Z, r ^ 1, then may be called as the secondary ones. Some of the series of these secondary 
poles may be truncated due to the particular sequence of the residues: 1, . . . , n, —n, . . . , —1, 0. □ 
The "method of steps" described above does not work when \q\ = 1. 

When |g| <C 1 it is possible to construct the solution of (JTJP by expansion of f(z,q) into a 
perturbation series over q. Plugging f(z, q) = J2uLq ^i{z)q 2 \ /(0, q) = 0, into QHD and solving the 
resulting differential equations for hi(z), one obtains the function /. For the first non- vanishing 
correction we have (/i = 1): 

f(z, 0) = h (z) = tanh z, 

h[{z) + 2h l (z)tanhz + l = 0, hi(0) = 0, 

2 z + ^smh2z 

j (z, q) = tanh z — q 2 h . . . . 

2 cosh z 

As it was shown above, the properties of the self-similar potentials differ for q being a root of 
unity of odd and even degree. We study below the first nontrivial cases, q 3 = 1 and q A = 1, in 
detail. 



4. Exact Solution for q 



3 _ 



Let us first derive a general solution of the dressing chain at iV = 3-periodic closure, which 
evidently contains as a particular case the g 3 = 1 self-similar solution. System (Q) consisting of 
three equations has two invariants of motion generated by TrA = —XIq — I\\ 

J = 2(/ (x) + /i(z)+/ 2 (aO), (21) 

h = kEf!-E**fj-*$- ( 22 ) 

Deriving explicitly the matrix A one finds (3{x) = ^ofi + /0/2 + Ai — A = 7(2;) — A. Substituting 
/3(x) in this form into (0) and setting the spectral parameter A equal to 7(2) we obtain: 

-i(7') 2 = (7-^o)(7-^i)(7-^)- (23) 
This is the differential equation for a Weierstrass p-function | 15| : 

(p') 2 = Ap 3 -g 2 p-g 3 , (24) 

p(x + 2u) = p(x + 2d) = p(x), p(x) = ^ + ||x 2 + ||x 4 + . . . , 
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where uj and uj' denote complex semiperiods. 

The expression for the potential uq is obtained by setting equal to zero an expression in front 
of A 2 in (jig): 

2 2 

u (x) = -2j(x) + J] Ei = 2p{x + x ) + |E ^. (25) 

i=0 i=0 

In the self-similar case one has Iq = independently on the value of a = /o(0). Without loss 
of generality we can set ciq = 0, which fixes I\ = 0, i.e. TrA = 0. Constants Ei are just equal to 
Aj: Ei = Aj. Setting Xj = -T^-j-g 2 - 7 , and choosing fi = q — 1 we get real Uq(x) for real x: 



e'2 



u (x) = fl - + A Q = 2p(x + lu 2 ), uo(0) = -2Ai = 2, (26) 

where is the equianharmonic Weierstrass function and uo 2 is the corresponding real semiperiod, 
oj 2 = uj + uj' . This p-function is defined by the conditions g 2 = 0, g 3 = 4 and it has definite ratio 
of semiperiods leading to very simple transformation rules under the scaling of its argument by q: 

to' = e 2m/3 uj, p(qx) = qp(x). (27) 

The fi(x) can be determined by solving either the differential equation (p6|) , or the system of 
three algebraic equations fl2lD, (p^), and (|25|): 

/t (x) = - ? ^ + ^) (28) 

where we put uji = qu 2 , u 3 = q 2 u 2 , LJ i+3 = LUi. 

The potential u (x) is periodic with minima at the points x\ = 2lu 2 and double poles at the 
points a;; = (2/+l)a;2, / G Z. The other two potentials are: U\{x) = 2p(x+uj 3 ), u 2 (x) = 2p(x+LUi). 
Note that all fa are complex and only uq is real. In the Appendix we present exact solution of the 
spectral problem for uq defined by the zero boundary conditions at singular points. 

5. The Case q 4 = 1 

When q 4 = 1, q = ±2, the general solution of the equation (|TJ) looks as follows: 

2/ W=ffW + ^-^|, (29) 

where g(z) is an arbitrary function subjected to simple constraints: 

g(qz)g{z) = fiq, g{Q) = (1 + q)f(0). (30) 



The conditions (p0|) for g(z) are compatible when /(0) = ±y/jlq/(l + q), i.e. a solution exists only 
for a definite value of /(0). In particular, there is no function g(z) such that g(0) = when \i ^ 0. 
The family of functions g(z) which satisfy the functional equation (|30|) is rather rich. We 



consider here two of the possible representations. The first one is: 

g(z) = ±^exp{F(z 2 )}, (31) 



S 



where F(y) is an arbitrary odd function of its argument (y = z 2 ). Corresponding potentials in 
the Schrodinger equation ([!]) are real when JJIq and F are real: 

«M = ™ + V + + f **2F<y) =F (32, 



When -F(y) — > +00 at 2; — > ±00 and /x 7^ 0, potentials d32|) grow exponentially and, hence, have a 
purely discrete spectrum. In the limiting case \i — > the choice oc ?/ yields the potential of 
a harmonic oscillator. 

The second solution is related to the (pseudo)lemniscatic Weierstrass function [T5|] which cor- 
responds to the invariants #3 = 0, #2 = ±1- Explicitly, 

p /2 = 4p 3 ±p, J = tuj, p(±iz) = -p(z), (34) 

u(z) = 2p(z + Zq), z = ^(. 

Taking z, zq to be real we obtain again the simplest Lame equation with the real potential. The 
corresponding spectral problem is exactly solvable (see Appendix). 

Note that for a general p-function and arbitrary q one has p(qz; qui, qui') = q~ 2 p(z; u>, ui'). The 
p-function constructed by the periods qui, qui' is equal to that built up by if the lattices of 
periods T = {u,u'} and V = {qui, qui'} are equivalent. But it is known that for two-dimensional 
lattices the only possible groups of rotations preserving the lattice are C n with n = 1, 2, 3, 4, 6. So, 
the transformation properties (p^) and (|34]) of the Weierstrass jp-function obtained for the cases 
q 3 = 1 and g 4 = 1 can take place in such a form only for one more non-trivial case when g 6 = 1. 
It is possible to extend the analysis to the higher roots of unity when the potentials are given by 
(O), using the appropriate transformation properties of the Riemann B-function. 



6. Realization of the q-Oscillator Algebra 

Group-theoretical content of the self-similar potentials discussed in the previous sections is de- 
scribed by the g-deformed Heisenberg-Weyl, or g-oscillator algebra ||: 

a~a + — q 2 a + a~ = /i, [a 1 * 1 ,/!] = 0, (35) 



where a are g-analogs of the creation and annihilation operators. The explicit realization of (j35|) 
utilizes the scaling operator T q , 

TJ(z) = f(qz), T q ± = -^T q , (36) 



T, ' f ' ' f ' — 1 rri rji -1 

qJ-r -Lqri q Q~ 1 > 1 ' 



and looks as follows 



a+ = ^ { ~Tz + q))Tq ' a ~ = ^ 1{ Tz + (37) 
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where f(z, g) is a solution for the equation fll4|) at x — z/q. For simplicity we keep in this section 
[i > 0. 

When g 2 — > 1 the relations ([$5|) are formally reduced to the standard bosonic oscillator com- 
mutation relations. As it was shown above, our system has a well-defined limit q — > 1 for arbitrary 
initial condition (|15|) ; however, the q — > —1 limit exists only when f(0,q) = 0. Although in both 
cases we get a harmonic oscillator potential, the second case corresponds to the non-standard 
realization of the Heisenberg-Weyl algebra: 

b + = ±i(--^ + ^z)P, b- = TiP(^ + \liz), [b-,b + }=fi, (38) 

where P is the parity operator, Pf(z) = f(—z), P 2 = 1. This observation displays the abil- 
ity of the g-deformation procedure to connect continuously different realizations of the original 
undeformed algebra. 

Consider the limit q — > 0. Then, one has formally a~a + = //. This algebra admits large variety 
of realizations, e.g. a~ may be a I x m rectangular matrix. In our case this limit is singular; 
however, the relation a~a + — q 2 a + a~ = fi still may be meaningful. In particular, if we substitute 
z = qx into (j37[) , then lim g ^ q 2 a + a~ = —d 2 /dx 2 ^ and lim^o a~a + — — d 2 /dx 2 (if z is kept 
finite, then g — > corresponds to the limit x — > oo which is not analyzed in this paper). Although 
in this case the self-similar potential is reduced to the known exactly-solvable one, we exclude the 
point q = from the further discussion because operators (|3~7D are not well-suited for the algebraic 
treatment of corresponding systems. 

In the g-deformed case the notion of number operator is not universal. In order to show this, 
let us consider the following combination of the ladder operators a ± : 

L = a + a- - v (39) 
= -^ + f(z)-f'(z)-u, " 



dz 2 1 — q 2 ' 

where we assume that q 2 ^ 1. Operator L satisfies the relations 

La ± = q ±2 a ± L } (40) 

which form a g-analog of the spectrum generating algebra of a harmonic oscillator problem, 
[N, b^] = ±b ± , N = b + b~. However, equations (|40|) do not imply the existence of a well-defined 
number operator iV with the properties [N, a^] = ±a ± . Indeed, the algebra (|3~5D, ([40]) has finite- 
dimensional representations for which any power of does not vanish (see below). The latter 
means that the spectrum of N would be unbounded neither from below nor from above, which 
in turn would contradict the finite-dimensionality of these cyclic representations. Note that usu- 



ally the relations (40) are used for the definition of the formal operator L itself. Sometimes the 
g-oscillator algebra is defined as a set of identities involving a^, //, L, and the inverse L~ l as well. 
We do not assume invertability of L, and, in a sense, utilize the minimal version of the g-oscillator 
algebra (|35|), when the operator L is expressed through a as it is given in d39|) . 



Let us discuss briefly the representation theory of (0), (|40| ) paying most attention to the 

(r) 

unitarity and finite-dimensionality of modules. Denote by ip\ the eigenstates of the abstract 
operator L, 

L^ = \4\ r = l,...,d, (41) 
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where A is some (complex) number. On the algebraic level nothing can be said about the degree 
of degeneracy d. For our explicit realization (|4lD is a standard Schrodinger equation which has 

(r) 

two independent solutions, i.e. d = 2. We suppose that the states tp x are uniquely fixed by some 
(boundary) conditions for all A. Relations ( f40|) allow to write 

« + # = « r + s (Ag 2 ) ^% = a" (A) ^« (42) 

Substituting this into (|35|) and using definition ([39|) we deduce 

«±(A) = a^(A)M r ± s 1 (A), a+(A)a~(A) = v + A, (43) 

where M rs is some non-degenerate matrix. The exact form of M rs (A) and a^(A) depends on the 
definition of the basis vectors ij) x r ' and can not be found algebraically. 

Proposition 2. There are only five types of essentilaly different from each other unitary modules 
of the algebra © (q ^ 0) arising at: 1. A < 0, < q 2 < 1; 2. A > 0, < q 2 < 1; 3. A = 0, 
q ^ ±1; 4. A ^ 0, -1 < g 2 < 0; 5. A ^ 0, g 2 = -1. 

Proof. We call the module unitary when the operators are hermitian conjugates of each 
other in the corresponding basis of states. The hermiticity conditions are ensured by A, q 2 being 
real and the choice a+(A) = (a~(A))* = \Jv + A, = M~ l . Let us consider first the interval 
< q 2 < 1. For negative eigenvalues A the representation should be truncated from below; 
otherwise v + A will start to be negative. So, we get the highest-weight infinite-dimensional 
representation: 

m = -»q 2j \j), 



a + \j) = f(l-gW)|j + l>, J = 0,1,... 
a-\l) = yjv{l-q^)\l-l), Z = l,2,... 

o"|0) = 0, (44) 

which is generated by the vacuum state |0) corresponding to the particularly chosen ipll. Note 
that we absorbed the matrix M rs ( |4*3| ) into the definition of states \j). In our case vacuum state 
is unique because it is defined by the first-order differential equation. The eigenvalues Xj = —vq 2 i 
provide the physical discrete spectrum of the infinite-soliton system described by our self-similar 
potential at real z, < q 2 < 1, \i > @, |3|]. The module ( |44"D is well-defined when q 2 — > 1. 
For positive A we get a non-highest-weight infinite-dimensional unitary module 

L\j)x = M 23 \j)x, jeZ 

a + \j)x = yJu + \q*U+V\j + l) Xl 



a~\j)x = sjv + \q 2 i\j-l)x. (45) 

Here we assume that the module is generated by some particularly chosen eigenstate |0)a of the 
Hamiltonian L with eigenvalue A, and we absorbed again the matrix M rs into the definition 
of states \j)x- Note that in this case a + lowers the energy. For our specific realization of the 
g-oscillator algebra wave functions form a particular subset of the scattering states of the 
infinite-solitonic potential. The limit q 2 — > 1 is not defined for (^5|). 
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A peculiar situation is described by the subspace of states corresponding to zero eigenvalue of 
L. In this case a + and a~ represent the integrals of motion, [a + ,a~] = [L, a ] = 0, and they can 
be diagonalized simultaneously with L: 

L^to = 0, a^ r) = e ±l6r ij (r) . (46) 

This is a d- dimensional (for any real q 2 , q 7^ ±1) reducible representations. The irreducible one- 
dimensional components are the simplest cyclic representations for which the number operator 
can not be defined. In our case d = 2 and in principle all angles 8 r may be found explicitly. The 
wave functions ip^ provide the simplest coherent states of the g-oscillator algebra. 

Finally, when — 1 < q 2 < we have the following possibilities. At |g| < 1 the infinite- 
dimensional representation with the defining relations as in ((44]) arises, but now the successive 
eigenvalues of L have different signs. At q = ±2 we have the following finite-dimensional repre- 
sentation: 

L\Q) X = A|0) A , L|1) a = -A|1> a , 

a+|0) A = y/17=\\l) x , a+|l) A = v / ^+A|0) A , 

a-\0) x = V7+X\0)x, ar\l) x = \/^A|0) A , (47) 

where A belongs to the interval [— v, v\. In this case the g-oscillator algebra is reduced to the 
following superalgebra: 

{a",a + } = /i, {a ± ,a ± } = /i±, [fj,, a ± ] = [/i±, a T ] = [/i±, y\ = 0. (48) 

The operators /x± may be set equal to zero, which corresponds to the boundary values of the 
parameter A in (f4TD, A = ±v, and this gives the standard fermionic oscillator algebra. □ 

Discussion of the physical significance of the states fl46f) , and analysis of the structure of 
modules provided by the self-similar potentials at — 1 < q 2 < lies beyond the scope of the 
present work. 

It is easy to see that the representations can be finite-dimensional iff X(q n — 1) = 0, q 7^ ±1. 
Indeed, the ordinary bosonic oscillator algebra arising at q — ±1 does not have finite-dimensional 
representations; the possibility q = was already excluded due to the exoticity. The case A = 
was considered above and it corresponds to the <i-dimensional, in general reducible, representation. 
At A 7^ all modules are constructed from series of states ^AgL for some range of j. Let q be 

not a root of unity. Then the set r \, } can be finite-dimensional iff it is truncated from below 

and above due to the zeros of af s (X). The equation a'lp^ = holds at A = —v and a + ^ A ^ = 
at A' = —vq~ 2 . Evidently, these states belong to two different irreducible infinite-dimensional 
highest-weight representations. We thus conclude that q n = 1, q 2 7^ 1 is a necessary condition for 
the existence of finite-dimensional representations at A 7^ 0. 

Finite-dimensionality of the modules at q n = 1 (n is the lowest number satisfying this identity) 
emerges due to the existence of the non-trivial central elements of the algebra: 

= (±ta ± )^ n \ [C ± , a T ] = [C*, L] = 0, (49) 

where [n] = n for odd n and \n\ = n/2 for even n. Due to the definition of L the states 
■0Ag 2 J) 3 — 0, . . . , \n\ form the \n\ -dimensional irreducible module which can not be unitary at 
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q 2 7^ —1. The non-zero C ± correspond to the cyclic representations ("cyclic" oscillators) for 
which the number operator does not exist. By purely algebraic means one finds that C satisfy 
the following polynomial algebra: 

C ± C ^ = V M -(-^H, (50) 



which is the particular case of algebras, considered in [II], |T2| from a different point of view. A 



subset of states, for which the right hand side of ([50]) vanishes, form the highest weight finite- 
dimensional representation. Note that the operator L^ n l enters (^) with different signs for n = 4k 
and n = 4k + 2. 

Returning to our explicit realization of the g-oscillator algebra we note that for odd n the C^ 1 
are ordinary differential operators of the n-th order whose commutativity with the Hamiltonian 
defines the particular cases of hyperelliptic potentials (|T2|) . In general the operators are not 
equal - this depends on the choice of the integrals of motion. E.g., the q 3 = 1 system fl28|) gives 

c+ = c~ = c, 

For the spectral problem considered in the Appendix, the L 3 + 1 operator is self-adjoint, but it is 
easy to see that its square root C does not preserve imposed boundary conditions. One cannot, 
however, exclude the possibility that for some of the spectral problems the operators C ± represent 
real physical observables. 
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Appendix 

Though the spectrum of a particle described by the Schrodinger equation 

- ij>" + 2p{x)^j = \if>, ^(0) = V^(2cu 2 ) = 0, (52) 
is known (see, for example 0]), we find it to be useful to show here how it can be derived from 



the chain (|2]). This derivation serves as an application of the factorization method to the spectral 
problems associated with singular potentials given by the (hyper)elliptic functions which appear 
after the periodic closure of the dressing chain. 

It was found that F(x) = ±w 2 ^ x ^ is the solution of the Riccati equation (|TTD , related to 
the Schrodinger equation (j52|). Taking x = to be a singular point of (3{x) and plugging in 
(3{x) = —p(x) — A, we obtain: 

_ F(x) _ t - ±w ± ( 53 ) 
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where w 2 = (trA) 2 — AdetA = — 4(A 3 + 1) for the q 3 = 1 system. Integrating (|53"D, one obtains the 
general solution of (0): 



ip(x) 



( r w — p'\ f r w + p' 



p + Asin (vWTlJ* + * I • (54) 



The quantization condition ?/;(0) = yields 5 = 0. When x — > 0, -y/p + A ~ x 1 , but sin (J 3 p ^ +A ; ■ ■ 
a; 3 , so that eigenfunction ip(x) ~ x 2 as a; — > 0. The condition ?/>(2u; 2 ) = yields the transcendental 
equation for A which determines the spectrum: 



v / A T TT / . " - = |7m. (55) 



The A — > oo asymptotics is: 



\M- 7r n , . _ 1 f + °° dx ln/1 1^ 191 

This is the spectrum of a particle in the infinitely deep well, as might be expected. For the first 
three eigenvalues numerical solution gives: \^ = 2.34, A^ 2 -* = 5.39, A® = 9.18. For any A^ one 
can find that if) ~ (x — lio-i) 2 when x — > 2co>2. 

This method of derivation of spectrum is not restricted to the equianharmonic Weierstrass 
function (self-similar case); the straightforward generalization of fl55"D yields the spectrum for 
arbitrary jp-function (El)f| 



n - s / ^t-t = \^ n( A - *y = a 8 - i^a - 93 ). (56) 

i=0 P\ x ) + A i=0 

The procedure presented above is not applicable in general for the g 4 = 1 potentials, because 
in this case the matrix element a±2 of the product of four Darboux transformations vanishes 
identically on the solutions of (|2]). However, the spectra of the particular cases corresponding to 
the (pseudo)lemniscatic Weierstrass function can be found from the formula ([56]). 

The obtained purely discrete spectrum (|5^) originates from the requirement of the finiteness 
of wave functions in the singular points of the potential. It is possible to remove singularities from 
the real axis by shifting the argument of the p-function along the imaginary axis. This leads in 
general to a complex potential. In particular, one cannot get a real potential without singularities 
in the case of the equianharmonic and pseudolemniscatic p-functions (when two of the roots Ei 
of the polynomial 4x 3 — g2X — g^ are complex). For the lemniscatic Weierstrass function the shift 
of the coordinate z — > z + u' leads to the real periodic bounded potential (— ~ < p < 0) with 
one finite gap in the spectrum. So, for the latter case one has two physically different self-adjoint 
spectral problems. 



8 The choice of constants Ei is restricted to the case when the corresponding potential is real and has singularities 
over the real axis. 
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